Mon. Not. R. Astron. Soc. 000, [T]j9] (0000) Printed 20 August 2012 (MN MfeX style file v2.2) 



Transit Timing Observations from Kepler: VII. 
Confirmation of 27 planets in 13 multiplanet systems via 
Transit Timing Variations and orbital stability 

Jason H. Steffen 1 , Daniel C. Fabry cky 2 ' 3 , Eric Agol 4 , Eric B. Ford 5 , 
Robert C. Morehead 5,6 , William D. Cochran 7 , Jack J. Lissauer 8 , 
Elisabeth R. Adams 9 , William J. Borucki 8 , Steve Bryson 8 , Douglas A. Caldwell 10 , 
Andrea Dupree 9 , Jon M. Jenkins 8 , Paul Robertson 7 , Jason F. Rowe 8 , 
Shawn Seader 8 , Susan Thompson 8 , and Joseph D. Twicken 8 

1 Fermilab Center for Particle Astrophysics, P.O. Box 500, MS 127, Batavia, IL 60510 

2 U CO /Lick Observatory, University of California, Santa Cruz, CA 95064, USA 

3 Hubble Fellow 

4 Department of Astronomy, University of Washington, Seattle, WA 90195 

5 Astronomy Department, University of Florida, 211 Bryant Space Sciences Center, Gainesville, FL 32111, USA 

6 National Science Foundation Graduate Research Fellow 

7 McDonald Observatory, The University of Texas, Austin TX 78730, USA 
8 NASA Ames Research Center, Moffett Field, CA, 94035, USA 

9 Harvard- Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA 
10 SETI Institute, Mountain View, CA, 94043, USA 



20 August 2012 



ABSTRACT 

We confirm 27 planets in 13 planetary systems by showing the existence of statistically 
significant anti-correlated transit timing variations (TTVs), which demonstrates that 
the planet candidates are in the same system, and long-term dynamical stability, which 
places limits on the masses of the candidates — showing that they are planetary. All 
of these newly confirmed planetary systems have orbital periods that place them near 
first-order mean motion resonances (MMRs), including 6 systems near the 2:1 MMR, 
5 near 3:2, and one each near 4:3, 5:4, and 6:5. In addition, several unconfirmed planet 
candidates exist in some systems (that cannot be confirmed with this method at this 
time). A few of these candidates would also be near first order MMRs with either 
the confirmed planets or with other candidates. One system of particular interest, 
Kepler-56 (KOI-1241), is a pair of planets orbiting a 12th magnitude, giant star with 
radius over three times that of the Sun and effective temperature of 4900 K — among 
the largest stars known to host a transiting exoplanetary system. 

Key words: celestial mechanics; methods: data analysis; techniques: photometric 



1 INTRODUCTION 

NASA's Kepler mission continues to identify many candi- 
date transiting exoplanet systems, which now number nearly 
2000 (|Borucki et al.||2011| |Batalha et~aT1|2012|. The pro- 



cess of confirming or validating these planet candidates as 
real exoplanets often requires a significant amount of anal- 
ysis and effort; and consequently the number of confirmed 
planets lags far behind the new candidate discoveries. In an 
effort to ameliorate this situation we published a new ap- 
proach to planet confirmation that requires a less detailed 



analysis (|Fabrycky et aT1|2012a| |Ford et al.||2012a| |Steffen 
et al.|2012a ) and can generally be accomplished in much less 



time. This approach essentially relies on demonstrating that 
two transiting candidates are in the same system and that 
their masses are planetary. Both aspects of this confirma- 
tion method rely on dynamical interactions. First, planets 
are shown to be in the same system by looking for anti- 
correlated transit timing variations (TTVs) that arise from 
short-term changes in obital period due to planet-planet in- 
teractions within the system ( Agol et al.|2005 Holman and 
Murray 2005). Second, the candidates are shown to have 
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planetary mass by requiring that the system be dynamically 
stable. 

This method is a particularly useful means of identi- 
fying true planetary systems from among the false positive 
systems — allowing ground-based follow-up resources to be 
devoted to studying systems that are both real and dynam- 
ically interesting (as manifest by their significant TTV sig- 
nal). While detailed TTV analyses have been instrumen- 
tal in confirming several multiplanet systems (e.g., Kepler 9 
(?), Kepler 11 ( |Lissauer et al.||2011a ), Kepler 18 (Cochran 



Table 1. Table of stellar properties for the planetary systems 
(taken from|Batalha et aT] (|2012|)). 



|et al. 2011), and Kepler 36 (Carter et al. 2012)), 21 plan- 



ets in 10 multiplanet systems have been confirmed by these 
new methods that we apply again here (Kepler systems 23 
through 32,|Fabrycky et al.|fl2012a| ); |Ford et al.| ( |2012a| ); [SteT 
fen et al. (2012a)). These 10 confirmed multiplanet systems 



include several that have particularly interesting properties 
and are likely to be the subject of future investigations. For 
example: all of the systems are near mean-motion resonance 
(MMR), Kepler-25 (KOI-244, [Steffen et al.| fl2012a|)) is rel- 
atively bright at Kp - 10, Kepler-30 (KOI-806, |Fabrycky 
et al. (2012a|)) has TTVs that deviate from a constant pe- 



riod by nearly a day over the course of the Kepler data, 
and many (Kepler 23, 24, 26, 31, and 32) have additional 
planet candidates that may be confirmed or validated. As 
more data from Kepler are gathered, new systems that can 
be confirmed by the means outlined above can be found. In 
this paper we apply the same methods used in (|Fabrycky| 
et al.|2012a]|Ford et al. 1 20 12a| [Steffen et al.|2012a| ) to 2 ad- 
ditional quarters of Kepler data (all data through Q8) and 
confirm 27 new planets in 13 systems. 

This paper is organized as follows. In ^2] we summarize 
the stellar properties for the host stars in these systems and 
the basic orbital and physical characteristics of the planets. 
The dynamical confirmation of the relevant 27 planets is pre- 
sented in ^3] We make concluding remarks in Q Data and 
results from some of the analyses performed in this paper 
are shown in the Appendix. 



2 STELLAR AND PLANET PROPERTIES 

We take the stellar properties from the tables given in 
|Batalha et al. (2012). The procedure used for determining 



the stellar parameters (effective temperature Teff, surface 
gravity \ogg and stellar radius R*) for our systems was dis- 
cussed in detail in Section 5.2 of that paper. In most cases, 
Teff and log g that were determined photometrically and re- 
ported in the Kepler Input Catalog were used as a start- 
ing point. These values were then corrected by matching to 
the Yonsei-Yale stellar evolution mod els (|Demarque et al~ 
2004j) (cf. Figure 2 of [Batalha et ah] fl2012|)). Kepler Fol- 



lowup Observing Program reconnaissance spectra have also 
been obtained for most of these KOIs. In general there is 
reasonable agreement between the stellar parameters in Ta- 
ble 1 and those obtained through the reconnaissance spec- 
troscopy. In a few cases of disagreement, the discrepancy can 
be attributed to the low S/N of the reconnaissance spectra, 
thus we rely on the published parameters. The basic planet 
size and orbital properties for the systems under considera- 
tion are also taken from (Batal ha et al.|2012 ) and are given 
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in Table [^f] We note that the giant star KOI- 1241 along with 
KOI-262 are quite bright and therefore amenable to in-depth 
asteroseismic study. We expect that the stellar properties, 
particularly for the giant, may be revised somewhat from the 
values presented here and in |Batalha et al.| ( |2012| once such 
studies have been completed. In neither case is the planet 
interpretation of the transiting candidates in jeopardy. 



3 CONFIRMATION VIA TTVS AND ORBITAL 
STABILITY 

3.1 Anticorrelated TTVs 



As mentioned above and discussed in some detail in (Fab- 
rycky et al.]|2012a| |Ford et al.||2012a| |Steffen et aT1|2012a| ), 
the presence of anticorrelated TTV signatures among planet 
candidates on a single target is a strong indicator that the 
objects are in the same system. Multiple methods can be 
used to establish these correlations — the approaches that 
we have used in the past include using Gaussian Process 
(GP) modeling (Ford et al. |2012a"| , Fourier analysis meth- 
ods ( |Steffen et al.||2012a| ), and the use of simple dynamical 
modeling of the potential planetary system ( Fabrycky et al. 



2012a). Each of these approaches relies on differing amounts 
of physical underpinning (dynamical modeling applying the 
most physics while GP the least). They also use different 
fundamental assumptions regarding the specific form of the 
TTV signal (GP being the most general of the three while 
dynamical modeling is the least). The Fourier method lies 
between the other two methods both in terms of the as- 
sumed properties of the TTV signal and its physical under- 
pinnings. Consequently, these different methods are quite 
complementary — all are suitable for confirming some sys- 
tems but none are suitable for confirming all systems that 
lie within the reach of this overall approach. 



The GP method of |Ford et al.|fl2012a| models the TTV 
signal as a Gaussian process, deriving a generalized func- 
tional form for each object in a given pair. From these GPs 

1 We note that because these planets show significant TTVs, the 
estimates for the transit depth and duration will be off slightly 
since the Kepler pipeline assumes a linear ephemeris. 
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Table 2. Table of planet orbital properties for the planetary systems. 
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the correlation between TTV signals of the two planet can- 
didates is calculated. Finally, we estimate the probability 
that the observed anticorrelation is due to timing noise with 
a Monte Carlo analysis. For this analysis we scramble the 
TTV residuals, remove any resulting linear trend, model the 
new residuals as a GP, and recalculate the correlation coef- 
ficient. We require, as a threshold to claim that the correla- 
tion is dynamical in nature, that less than 0.1% of the Monte 
Carlo realizations have a more significant correlation than 
the actual data. 



on a grid of periods: 

/i = Asin(^)+Bcos(^)+C. (1) 

where A, B, and C are model parameters, and Pi is the test 
timescal^] The fitted values for A and B are stored along 
with their measured uncertainties, a a and as, derived from 
the covariance matrix of the three-parameter fit. Next, a 



The Fourier method of |Steffen et al. (2012a) generates a 
periodogram of the TTV signal by fitting sinusoidal models 



2 This formula, originally given in |Steffen et al.| (|2012a >, had a 
typo that has been corrected here and in |Steffen et al.| ( 2012b| >. 
Specifically, P I Pi is replaced by 1/Pi in the argument of the si- 
nusoid functions. 
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quantity H 



terns (e.g., |Ford et al.| (|2012b|); |Steffen et al] {[2012b)). In 



A X A 2 B X B 2 



&A1CTA2 



CTB1CTB2 



(2) 



is calculated for each of the sampled periods where the "1" 
and "2" subscripts correspond to the two objects. Once these 
values have been calculated, the maximum value of S for a 
given candidate pair, S ma x, is recorded. Finally, a Monte 
Carlo analysis is done to identify the probability of ran- 
domly finding a S ma x as large as that observed. Again, our 
threshold for claiming evidence for dynamical interactions is 
0.1%. 



The dynamical modeling of |Fabrycky et aT (2012a) 



seeks to extract the dominant TTV signal from numeri- 
cal integration, and determine if this signal is present in 
the data. A numerical integration of the system is run for 
each pair of planets in Table [2] We assume circular, copla- 
nar orbits with the same orbital properties (period, phase) 
and physical properties (mass, derived from radiu s via the 
power- law M p — M®(R P /R { 



Lissauer 



as was used in 
et al. ( 2011b| ) as the planets we are modelling. From that 
simulation, the dominant frequency of the resulting TTV 
signal is found in a periodogram (Lomb 1976). We fit a si- 
nusoid of that frequency to the data, determining its ampli- 
tude. Again, a Monte Carlo study is done to see if a ran- 
dom realization can produce that large an amplitude at the 
theoretically-predicted frequency. One unique characteristic 
of the dynamical modeling, contrasted with the GP and 
Fourier methods, is that it does not require a significant 
TTV signal from both objects in the tested pair — and there- 
fore not requiring an anti-correlation in the timing residuals. 
Rather, we only use the transits of the perturbing planet to 
determine its period and phase, from which the perturbed 
planet's expected TTV signal follows. Thus, while it may 
lack generality concerning the types of TTV signals it can 
find, it compensates somewhat for this lack with its power 
to confirm the cohabitation of candidates in a system when 
one TTV signal is very weak — e.g., if a high-mass planet is 
being perturbed by a low-mass one. 

Each of these methods was applied to all of the systems 
being considered. The transit times for each of the analyzed 



candidates are produced in the manner described in (Ford 
et al.|2011| ). For each system, outlier data are removed by 
eliminating all transit times that are four times the Median 
Absolute Deviation from the best fitting linear ephemeris 
and eliminating all transit times where the timing uncer- 
tainty is more than twice the median timing uncertainty (as 
was done in |Ford et aT| ( |2012b| ) and |Steffen et al.| ( |2012b] )). 
(The timing residuals for the sytems are shown graphically 
in the Appendix. ) The results of the Monte Carlo studies, 
specifically the expected False Alarm Probabilities (FAP) 
for relevant cases with each method, are given in Table [3] 
In all cases at least one of the three methods generates a 
sufficiently low FAP to claim that the observed autocorrela- 
tion is not due to timing noise, but more likely two planet 
candidates in the same system. In many cases all methods 
can make this claim. 

In a few cases (planets in KOIs 886, 1336, 1529, and 
2086) the signal to noise ratio (SNR) of each individual tran- 
sit event is relatively small. Low SNR transits occationally 
result in poorly determined transit times and often these 
cases are ignored in general TTV studies of the Kepler sys- 



order to increase our confidence that the observed TTV sig- 
nal is real, we modelled the lightcurve using a straightfor- 
ward linear-plus-sinusoid ephemeris in which we allowed the 
times of transit to have a sinusoidal modulation. The transit 
duration, impact parameter, and depth for each transiting 
planet were fixed (since variations in these quantities have 
yet to be detected for any small exoplanet). In all cases we 
fit the transits of both planets simultaneously. The results 
of this expanded ephemeris yielded a significant reduction in 
the x 2 over a simple linear ephemeris model (the reduction 
in x 2 being 830 for KOI-886, 47 for KOI-1336, 130 for KOI- 
1529, 48 for KOI-2086 .01 and .02, and 166 for KOI-2086 .02 
and .03). These results demonstrate that the modulation in 
transit times, with its anticorrelation and timescale, is not 
due to poorly determined transit times and we proceed with 
our analysis in the same fashion as the other systems. 

3.2 Discussion of plausible false alarms 

The choice of an FAP < 10 -3 is the same as was done in 



Fabry cky et al.| ( |2Q12a| ), |Ford et al.] ( |2012a| ), and |Steffen 
et al. (2012a). We have not published an exhaustive search 



for TTV systems that can be confirmed with this method 
as some Kepler systems have warranted individual scrutiny. 
In addition, the vetting of all KOIs is still ongoing — even 
as the list of KOIs continues to grow — and therefore re- 
mains nonuniform. Nevertheless, we do want to be confi- 
dent that we are not misinterpreting these systems. Anti- 
correlated TTVs with a well-defined period and phase pro- 
vide a distinct functional form that yields high confidence in 
the pair being true interacting planets. However, autocorre- 
lations between parabolic TTV signatures imply a gradual 
change in the orbital periods, but are less powerful at elimi- 
nating alternative explanations. We look at two possible ex- 
planations for the observed anticorrelated TTVs and show 
that they are significantly less likely than the single sys- 
tem interpretation — justifying our claim that the observed 
TTVs are due to planetary interactions in a single system. 
First, we consider the possibility of a background star that 
hosts one of the planets. Second, we consider the case of a 
bound stellar binary with a planet orbiting each star where 
the light travel time across the system produces the TTV 
signal. 

Consider the case where both the target star and a back- 
ground star host a transiting planet; and that these planets 
are both being perturbed by an unseen planetary compan- 
ion. Occasionally a parabolic TTV signature from one star 
may be anticorrelated with a parabolic TTV signal from 
the other. To estimate a limit on the probability of this sce- 
nario (we do not attempt to give an exact estimate, only 
to constrain its value) we multiply: 1) the fraction of KOIs 
that show only a quadratic TTV signature, 2) the fraction 
of KOIs as near to first-order MMRs as those in our sample, 
3) the probability of a given target having a planet can- 
didate (a proxy for the fraction of background stars that 
would show a candidate), and 4) the number of candidate 
pairs that we investigated — then we divide by two since we 
require an "anti" correlation in the TTV signals. Recall that 
the resulting false alarm probability is (a limit on) the prob- 
ability of misinterpreting two distinct planetary systems as 
a single planetary system. 
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Using data from Ford et al. (2012b) we estimate the 



fraction of KOIs that show evidence for quadratic TTVs by 
selecting KOIs from ( |Ford et~aL 2012b| ) where a fit using 
a quadratic ephemeris has a x p-value less than 0.05. We 
also require that a cubic or quartic ephemeris fail the F-ratio 
test (with a p- value greater than 0.1) when compared to a 
quadratic ephemeris — we want only systems with quadratic 
trends as the others are not subject to this particular mis- 
interpretation of the TTV signal. We find that - 20% of 
Kepler planet candidates show some evidence for quadratic- 
only trends. Next, we note that the planet pairs we search 
for are within ~ 10% of a first-order MMR where we ex- 
pect TTV signatures to be manifest. The largest absolute 



may indicate the possibly of a false positive. As a check, we 
generated a synthetic population of planet canidates using 



a Monte Carlo code similar to iLissauer et al. (2011b ) and 



deviation from zero, using the £i statistic from (Lissauer 
et al.||2011b| ) and flFabrycky etaTl|2012b) , is 0.37. We esti- 



mate the probability of finding two randomly chosen planet 
candidates having £i < 0.4 using the distribution of KOI 



periods from (Bat alha et al.|2012 ). We select all planet can- 
didates with period ratios within a factor of three of the 
median period of all KOIs (essentially the peak of the dis- 
tribution which yields a conservative estimate) , calculate £i 
relative to the median period, and divide by the total num- 
ber of KOIs. The result is that no more than 15% of ran- 
domly chosen KOI pairs will be as close to a first order 
MMR as the candidates we study here. Finally, the number 
of KOIs seen by Kepler suggests that ~ 2% of Kepler tar- 
get stars will have a planet candidate. Thus, no more than 
(1/2) x (0.2) 2 x 0.15 x 0.02 ~ 1 x 10" 4 ) of randomly-selected 
pairs of candidates would show parabolic ant icorrelat ions. 
When identifying systems to study in this paper, we search 
~ 1000 pairs which yields a naive expectation of 0.1 candi- 
date pairs that show parabolic, anticorrelated TTVs. Com- 
paring this number to the ~ 10 previously confirmed systems 
that have have anticorrelated TTVs yields a less than 1% 
chance of misinterpreting planets orbiting two different stars 
as a system of planets orbiting a single star. 

Still, much more can be brought to bear to reduce 
this false alarm probability still further (not including con- 
straints on the population of suitable background stars in 
the galaxy). We look at two tests specifically. First, and most 
importantly, we examined the displacement of the centroid 
of the target during transit. If a pair of planets were orbiting 
different stars, then the centroid of the target would displace 
in different directions during the transits of each planet. In 
all cases the displacement of the centroid is consistent with 
the planet pairs orbiting the same star and that the star is 
indeed the Kepler target. 

For the second test, we calculated the quantity £ = 
(Di/ 'Z^Xft/ Pi) 1 ^ 3 for each planet pair where D\ and D2 
are the transit durations and Pi and P2 are the orbital pe- 
riods for the inner and outer planets respectively. The value 
of £ for each pair should be of order unity if the planets are 
in the same system ( Steffen et al.|2010 ). Indeed, the bulk of 
cases presented here are within the 98% confidence interval 
of the distribution for observed £ values for all of the Kepler 
multiple systems. The only exception is anomalous values in 
the Kepler-51 (KOI-620) system (2.57 and 0.38 for the 01 
- .03 and .03 - .02 pairs respectively). In fact, the £ value 
for the .03 - .02 pair is the lowest observed to date. In this 
system one planet, Kepler-51c (KOI-620. 03), appears to be 
on a nearly grazing orbit (see Batalha et al. (2012)) — which 



Fang and Margot (2012), drawing 2 to 10 planets per system 



uniformly over the Kepler candidate period distribution and 
assuming circular orbits. Radii were also drawn from the em- 
pirical Kepler distribution. The mutual inclinations of the 
planets with respect to an arbitrary reference plane were 
Rayleigh distributed wit h a mean of 1.52 degrees , equivalent 
to the best fit regions oflFabrycky et al.l d2012bk. After 10 5 



realizations, £ was computed for the synthetic candidates 
with multiple transiting planets having total SNR > 7.1. 
A two sample KolmogorovSmirnov test gives a P-value of 
0.14 indicating the null- hypothesis that the synthetic and 
observed distributions come from the same parent distribu- 
tion can not be rejected at high significance. From the syn- 
thetic population we calculate the probability of £ < 0.38 or 
£ > 2.57 to be 0.016. Given that there are 365 systems with 
multiple transiting planets and 529 corresponding £ pairs, a 
true planet system showing such extreme values is not un- 
expected. The £ values for the different planet pairs are also 
given in [2] 

The second possible misinterpretation of an anticor- 
related TTV signal that we consider is planets orbiting 
two gravitationally-bound stars. Here, the light travel time 
across the system would produce the observed period varia- 
tions. In many cases this stellar binary interpretation is not 
consistent with the amplitude of the TTV signal — the im- 
plied orbital period is short (a few hundred days), meaning 
that the distance between the stars would be small, but the 
TTV signal is large, requiring a large separation between 
the stars. The worst-case scenario is when: 1) the binary 
has an orbital period equal to the duration of the data, 2) 
the phase of the orbit being such that the TTVs actually 
have the form of a cosine instead of a parabola, and 3) the 
second star being bright enough to allow the observed tran- 
sit depth (here we use a stellar twin) but essentially massless 
so that the light travel time is maximized — a somewhat un- 
physical worst-case scenario. None of the systems that can 
be modeled with sinusoids can be explained with the binary 
star model since their TTV amplitudes are too large. The 
amplitudes of all observed parabolic TTV signals are also 
too large by a factor of at least 7 with the exceptions of 
Kepler- 50 (which is very near the 6:5 MMR) and Kepler- 
60 (which has three planets instead of just two). Still, even 
these systems have amplitudes that are a factor of 2 and 1.5 
too large respectively for this worst-case scenario. 

In these remaining two systems (Kepler-50 and Kepler- 
60), all of the candidate pairs are very near first order 
MMRs — all are comfortably within 10 -3 of the resonance 
and have Ci < 0.035. Planets this close to MMR would 
naturally have TTV signals that develop over very long 
timescales (and would be seen only as parabolas over the 
course of our data). Using the same method outlined above 
we find that fewer than 2% of planet pairs are this close 
to resonance (having (1 < 0.04) — and for Kepler-60 one 
would need two such pairs. The probability of having a bi- 
nary star system with an appropriate period ratio — between 
500 and 1500 days — is approximately 0.05/2 = 0.025 (us- 
ing the analytic distribution given in Duquennoy and Mayor] 
( 1991 ) and dividing by two in order to have a binary star 



naturally produces such extreme values for £, however this instead of an isolated star). Also needed is that all plan- 
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ets must transit their respective hosts (with a probability 
perhaps as large as ~ 0.1 given that the binary star sys- 
tem which would likely induce correlations in the orbital 
planes of the planets), Finally, the orbital phase of the bi- 
nary must be such that it produces a parabola-resembling 
cosine signal 2 x 0.2). Thus, the expected number of 
planet pairs in an plausible, stellar binary system is no larger 
than 0.02 x 0.025 x 0.1 x 2 x 0.2 - 10" 5 — a result slightly 
smaller than the scenario using a background star. Again, 
this possibility is much less likely than that of a single system 
(nearly 1000 times less likely for these two systems where 
this model is potentially viable). Nevertheless, the twin- 
star scenario given here has always been a known caveat 
for multi-transiting systems (especially for very widely sep- 
arated binaries) and should such a system ever be seen, it 
would be a very interesting one for study. 

3.3 Mass limits from orbital stability 

Once the fact that the different planet candidates orbit the 
same star has been established, we use dynamical stability 
to show that their masses must be planetary. We followed 



Table 4. Mass limits from the stability analysis. 



the method used and described in Fabrycky et al. (2012a). 
Briefly, we integrate a series of planetary systems with dif- 
ferent masses, starting from the value at which stability is 
analytically assured ( |Gladmanfl 993 ) . We integrate the sys- 
tem until it becomes unstable, produces close encounters, or 
survives without close encounters for 10 9 inner-planet orbits. 
We set the upper limit on the planet mass as the small- 
est mass that yields an unstable system. For each planet 
candidate examined with the stability test, this value is 
< 25Mj up iter, confirming them as planets (using the planet 
criterion suggested by (Schneider et al. 2011| ). Table[4]shows 
the overall results of the stability study. 



Ford et al. (2012a) discusses in detail the possible, 
though unlikely, scenario where the planets in those systems 
are both orbiting a background star instead of the target 
star. We mention here the implications for the planet na- 
ture of the orbiting objects. Dynamical stability yields con- 
straints on the masses of the planets compared to the total 
mass of the system. So, if a background star is more massive 
than the target, then the corresponding mass limits would 
grow in proportion to the difference. It is possible that such 
a scenario would push the mass limits beyond the planetary 
regime. For most systems the mass of the true host star 
would need to be larger by a factor of 10 or more (for most 
of the systems here that would imply > 8M© host). Such 
massive stars are quite rare and don't survive on the main se- 
quence for very long. For the less massive targets with large 
mass limits on the orbiting objects — Kepler-52 (KOI-775) in 
particular — a more massive background star, with a factor of 
two or more greater mass, could be an issue (though the or- 
bits would all need to be oriented in a physically unfavorable 
manner to yield consistent transit durations). Nevertheless, 
refined mass limits from the next section effectively elimi- 
nate the possibility that the objects orbiting Kepler-52 are 
not planetary. 

3.4 Analytic mass limits 



Planet 


KOI 


Max Mass (Mj) 


48b 


148. 01 


5.94 


48c 


148. 02 


11.61 


49b 


248. 01 


0.98 


49c 


248. 02 


0.72 


50b 


262. 01 


0.10 


50c 


262. 02 


0.11 


51b 


620. 01 


3.23 


51c 


620. 03 


2.60 


52b 


775. 02 


8.70 


52c 


775. 01 


10.41 


53b 


829. 01 


18.41 


53c 


829. 03 


15.74 


54b 


886. 01 


0.92 


54c 


886. 02 


0.37 


55b 


904. 02 


1.49 


55c 


904. 03 


1.11 


56b 


1241.02 


5.12 


56c 


1241.01 


12.18 


57b 


1270.01 


18.86 


57c 


1270.02 


6.95 


58b 


1336.01 


1.39 


58c 


1336.02 


2.19 


59b 


1529.02 


2.05 


59c 


1529.01 


1.37 


60b 


2086.01 


0.25 


60c 


2086.02 


0.56 


60d 


2086.03 


0.68 



A recent paper by Lithwick et al. (2012) gives analytic for- 



mulae that can be used to derive mass estimates or upper 



limits using the relative amplitudes and phases of the TTV 
signal. Those formulae apply once a reliable amplitude for a 
sinusoid-shaped TTV signal can be extracted. Here we as- 
sume that the free eccentricity is zero and apply the corre- 
sponding formula from ( Lithwick et al.|2012 ) to the systems 
where this analysis is appropriate (i.e., where the amplitude 
of the TTV signal can be extracted). 

The free eccentricity of a planet — the eccentricity it 
would have in isolation (as opposed to the forced eccen- 
tricity, which results from the presence of the perturbing 
planet) — will tend to boost the TTV signal on a nearby 
planet. Thus, if there is some free eccentricity, then the mass 
estimates that come from the above analysis should be in- 
terpreted as approximate mass upper limits. A constraint 
on the free eccentricity can be made from the phase of the 
TTV signal relative to the line of sight of the observer. In 
our sample of systems the uncertainty in the TTV phase is 
sufficiently large that we do not claim absence of free eccen- 
tricity in any system. Thus, we consider all estimates to be 
upper limits. Nevertheless, these upper limits are roughly 
one order of magnitude more stringent than what one ob- 
tains from dynamical stability. 

The results of this analysis are presented in Table [5] 
where we include the mass limits from orbital stability 
(copied from table [4] but now in Earth masses), the analytic 
mass estimate assuming zero free eccentricity, the formal un- 
certainty in this estimate, and the measured TTV phase and 
its uncertainty (note that we quote the best fitting value for 
the masses, so the actual mass limit at a specific confidence 
level must be calculated from this quantity and its given 
uncertainty). The systems analyzed in this manner include 
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Table 3. Results of the Monte Carlo studies of the TTV residuals. Only pairs with an FAP < 0.1 for at least one of the methods are 
shown. If no realizations of the Monte Carlo were more significant than the observed signal out of 10 4 realizations, we note the results as 
FAP < 10 — 3 . Here FAPi and FAP2 come from the dynamical modelling where the inner planet is being pertubed by the outer and vice 
versa. FAP gp comes from the Gaussian Process Monte Carlo simularion. And, S and FAP^ are for the Fourier analysis Monte Carlo. 
Confirmed planet pairs are noted. 



Kepler 


KOIi 


KOI 2 


FAPi 


FAP 2 


FAP GP 


"max 


FAPs 


Confirmed 


48 


148.01 


148.02 


0.0022 


< 10- 3 




19.6 


< 10- 3 


y 




148.01 


148.03 
















148.02 


148.03 














49 


248.03 


248.01 
















248.03 


248.02 
















248.03 


248.04 




0.0816 












248.01 


o a no 
z4o.Uz 


< 10 


< 10 — 3 




69.8 


< 10 


y 




248.01 


O A Q C\A 








7.1 








248.02 


248.04 














50 


262.01 


ncn no 
ZOZ.Uz 


< 10 


0.0120 


< 10 


34.8 


< 10 


y 


51 


620.01 


620.03 


0.0001 


0.0166 


0.0006 


79.6 


0.0072 


y 




620.01 


620.02 
















con no 


con no 

bzU.Uz 


0.0114 












52 


775.02 


775.01 


0.0003 


0.0476 




19.1 


0.0001 


y 




775.02 


775.03 
















775.01 


775.03 


0.0439 












53 


829.02 


829.01 








5.9 








829.02 


829.03 




0.0019 


0.0004 










829.01 


829.03 


0.0009 


0.0007 




10.9 


0.0072 


y 


54 


886.01 


886.02 


< 10- 3 


< 10- 3 


< 10- 3 


104.3 


< 10- 3 


y 




886.01 


886.03 
















oou.uz 


QQf. no 
oou.uo 








u.u 






00 




nn/i nzL 








7 7 








904.01 


904.05 
















904.01 


904.02 
















904.01 


904.03 




< 10- 3 












904.04 


904.05 
















904.04 


904.02 




0.0001 












904.04 


904.03 








6.9 








904.05 


904.02 








6.2 








904.05 


904.03 
















904.02 


904.03 


< 10- 3 


< 10- 3 




93.2 


< 10- 3 


y 


56 


1241.02 


1241.01 


< 10- 3 


0.0210 




11.5 


0.0016 


y 


57 


1270.01 


1270.02 




< 10- 3 




18.1 


< 10- 3 


y 


58 


1336.01 


1336.02 


0.0001 


0.0001 




13.6 


0.0006 


y 


59 


1529.02 


1529.01 




< 10- 3 






0.28 


y 


60 


2086.01 


2086.02 


0.0002 


0.0092 




10.0 


0.0049 


y 




2086.01 


2086.03 






0.0001 


7.8 




y 




2086.02 


2086.03 


0.0194 


0.0001 








y 



Kepler-48, Kepler-49, Kepler-52, Kepler-53, Kepler-57, and 
Kepler-58. For these systems, the location in time where 
the line of conjunction is along the line of sight is shown in 
the corresponding figures in the Appendix. Zero TTV phase 
(measured for the inner planet) occurs when the TTV signal 
crosses zero from above. A more detailed study of these sys- 
tems and those discussed in ( |Lithwick et al.|2012 ) will yield 
an interesting estimate for the number of near-resonant sys- 
tems with free eccentricity. At first glance, it appears that 



rougly half of the systems studied using this approach re- 
tain some free eccentricity (~ 6 out of 12, though more data 
will inevitably yield more a more reliable count). Such a re- 
sult will have implications for how the excess of planetary 
systems near, but perhaps not in, low-order MMRs is gen- 
erate d ((jLithwick and Wu||2012 |Terquem and Papaloizou 



2007l|Papaloizou|2011||Podlewska-Gaca et al.|2012| |Sandor 
and Kley|2006| |Batygin and Morbidelli||2012| )) 
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Table 5. Mass limits from orbital stability (now given in Earth masses) and from analytic formulae with TTV phase cf). The quoted 
numbers from the analytic results are the best fitting values, so the actual mass limits will be a combination of these numbers and their 
uncertainties. 



Kepler 


KOI 


Stability Mass Limit 


TTV Mass Limit 




<$> (deg) 


0> 






(M©) 


(M e ) 


(M ) 


(degrees) 


(degrees) 


48b 


148.01 


1960 


17.2 


3.9 


-7.9 


14 


48c 


148.02 


3830 


10.1 


3.5 


-156 


10 


50b 


248.01 


323 


7.6 


1.3 


47 


9.0 


50c 


248.02 


238 


7.1 


1.1 


-138 


9.3 


52b 


775.02 


2870 


89 


46 


-153 


9.6 


52c 


775.01 


3440 


37.4 


8.0 


39 


27 


53b 


829.01 


6080 


84 


27 


-42 


26 


53c 


829.03 


5190 


24 


12 


111 


40 


57b 


1270.01 


6220 


100 


15 


98 


190 


57c 


1270.02 


2290 


5.4 


3.7 


-84 


60 


58b 


1336.01 


677 


27.4 


8.1 


-85 


122 


58c 


1336.02 


452 


41 


12 


135 


17 



4 DISCUSSION AND CONCLUSION 

All of the systems that we confirm here and in similar previ- 
ous studies are near first-order MMR. In this case there are 
six pairs near the 2:1 MMR (Kepler-48, Kepler-51, Kepler- 
52, Kepler-53, Kepler-56, and Kepler-57), five pairs near the 
3:2 MMR (Kepler-49, Kepler-54, Kepler-55, Kepler-58, and 
Kepler-59), and one pair each in the 4:3 (Kepler-60 c and d ), 
5:4 (Kepler-60 b and c), and 6:5 (Kepler-50). Several systems 
show additional planet candidates that, due to their pres- 
ence in known multi-transiting systems, are very likely to 
be true planets as well ( |Lissauer et al.|2012| . If we consider 
all of the KOIs to be planets, then a few systems will have 
multiple pairs of planets near first order MMRs with some 
of these pairs being adjacent links in a near resonant chain 
including: Kepler-51 near a 3:2 - 2:1 chain (3:2:1), Kepler-53 
near a 2:1 - 2:1 chain (4:2:1), Kepler-55 with one pair near 
the 2:1 and a separate pair near the 3:2, and Kepler-60 near 
a 5:4 - 4:3 chain (20:15:12). The balance of the planet pairs 
in these systems, if near a resonance at all, must be near a 
resonance of higher order. 

These planetary systems, like the other systems that 
have been confirmed via TTV analyses, are of particular 
interest for long-term scrutiny with follow-up observations 
and dynamical studies. The detailed dynamics of plane- 
tary systems yield meaningful information about the evo- 
lutionary histories of the orbital architectures of the sys- 
tem. For example, the differences in the fraction of systems 
that show multiple candidates and detectabe TTVs among 
smaller planet candidates with few-day orbital periods com- 
pared with the larger, Jupiter-sized candidates in the same 
period range — these hot Jupiters have no compelling signs 
of current planet-planet dynamics — gives strong clues as to 
their likely distinct dynamical history compared to the bulk 
of the exoplanet population (Stef fen et al.||2012| ). 

These newly confirmed planetary systems continue to 
show the value that transit timing variations have and will 
have in transiting planet endeavors. As new planets are 
found, with increasingly smaller sizes and longer orbital pe- 
riods, TTVs provide the only viable means of determining 
planetary masses. TTVs, along with statistical validation via 



the planetary nature of the planet candidates that the Ke- 
ller spacecraft identifies (in part because the stars in the 
Kepler field are often too dim to make good Radial Velocity 
targets) . 

Two of these systems show a modest gap in the orbital 
period ratios of adjacent pairs of KOIs. If planetary systems 



tend to be dynamically packed (e.g. Barnes and Raymond 
fl2004|); [Lissauer| (|1995); |Laskar| (f2000|)) then there may be 



the techniques employed in (Torres et al. 2011 ) and (Lissauer 
et al.|2012 ), will likely be the primary means to demonstrate 



additional undetected planets in these gaps. In this set of 
systems, the most promising examples for future investi- 
gations include Kepler-48 which has a gap that is over a 
factor of 4 in period ratio and KOI-248 which has a gap 
that is nearly a factor of three in orbital period. A detailed 
TTV analysis of these two systems would provide useful 
constraints on or support for the hypothesis of dynamically 
packed systems. 

The planetary system orbiting the giant star Kepler-56 
is another interesting individual system — being an early dis- 
covery of a transiting multiplanet system hosted by a giant. 
Chromospheric emission from giant stars causes the star to 
be limb brightened in certain pass bands (e.g. the Calcium 
H and K lines) . Comparison of the transit duration in these 
bands with the duration from the broad band Kepler data 
may yield a direct measurement of the size of the chromo- 
sphere ( |Assef et al.|2009 ). 

Continued monitoring by the Kepler spacecraft and fu- 
ture study will likely yield more precise mass measurements 
of the planets in these systems and consequently provide 
constraints on their bulk composition. Such mass measure- 
ments will in turn improve our understanding of planet for- 
mation processes outside of our own solar system and are 
useful to place or solar system in the context of the general 
planet population. 
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APPENDIX A: TIMING RESIDUALS FOR 
EACH SYSTEM. 

In this appendix we show the data that were used in the 
analysis of the paper (Figures 



Al 



A8 A9 A10 All A12 and A13). The timing residuals 



A2 



A3 



A4 



A5 



A6 



A7 



(TTV signal) of the inner planet (s) are displaced vertically 
for convenience in viewing. We also show in Table |A1| the 
S max values, the TTV period associated with S ma x, and as- 
sociated fit parameters (from Equation [I]) for the pairs of 
planets that can be confirmed with the Fourier methods. 

This paper has been typeset from a TgX/ WF^i file prepared 
by the author. 
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Table Al. Results from the EE max calculations for pairs of planets that can be confirmed, or nearly so, with the Fourier method. 
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Figure A2. Deviations from a constant period for Kepler-49 
(KOI-248). The residuals for the outer planet have been displaced 
vertically for convenience in seeing the TTV signal. Vertical lines 
correspond to the times when the line of conjunctions of the two 
planets crosses the line of sight and are used to measure the TTV 
phase for the analytic mass estimates shown in Table [5] 
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Figure A3. Deviations from a constant period for Kepler-50 
(KOI-262). The residuals for the outer planet have been displaced 
vertically for convenience in seeing the TTV signal. 
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Figure A4. Deviations from a constant period for Kepler-51 
(KOI-620). The residuals for the outer planet have been displaced 
vertically for convenience in seeing the TTV signal. 
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Figure A7. Deviations from a constant period for Kepler-54 
(KO 1-886). The residuals for the outer planet have been displaced 
vertically for convenience in seeing the TTV signal. 
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Figure A5. Deviations from a constant period for Kepler-52 
(KOI-775). The residuals for the outer planet have been displaced 
vertically for convenience in seeing the TTV signal. Vertical lines 
correspond to the times when the line of conjunctions of the two 
planets crosses the line of sight and are used to measure the TTV 
phase for the analytic mass estimates shown in Table [5] 
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Figure A6. Deviations from a constant period for Kepler-53 
(KOI-829). The residuals for the outer planet have been displaced 
vertically for convenience in seeing the TTV signal. Vertical lines 
correspond to the times when the line of conjunctions of the two 
planets crosses the line of sight and are used to measure the TTV 
phase for the analytic mass estimates shown in Table [5] 
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Figure A8. Deviations from a constant period for Kepler-55 
(KOI-904). The residuals for the outer planet have been displaced 
vertically for convenience in seeing the TTV signal. 
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Figure A9. Deviations from a constant period for Kepler-56 
(KOI-1241). The residuals for the outer planet have been dis- 
placed vertically for convenience in seeing the TTV signal. 
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Figure A10. Deviations from a constant period for Kepler-57 
(KOI-1270). The residuals for the outer planet have been dis- 
placed vertically for convenience in seeing the TTV signal. Verti- 
cal lines correspond to the times when the line of conjunctions of 
the two planets crosses the line of sight and are used to measure 
the TTV phase for the analytic mass estimates shown in Table [5] 
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Figure A 13. Deviations from a constant period for Kepler-60 
(KOI-2086). The residuals for the outer planet have been dis- 
placed vertically for convenience in seeing the TTV signal. 
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Figure All. Deviations from a constant period for Kepler-58 
(KOI-1336). The residuals for the outer planet have been dis- 
placed vertically for convenience in seeing the TTV signal. Verti- 
cal lines correspond to the times when the line of conjunctions of 
the two planets crosses the line of sight and are used to measure 
the TTV phase for the analytic mass estimates shown in Table [5] 
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Figure A12. Deviations from a constant period for Kepler-59 
(KOI-1529). The residuals for the outer planet have been dis- 
placed vertically for convenience in seeing the TTV signal. 
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